Anomalous Conical Di-jet Correlations in pQCD vs AdS/CFT 
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We propose an identified heavy quark jet observable to discriminate between weakly coupled 
pQCD and strongly coupled AdS/CFT models of quark gluon plasma dynamics in ultra-relativistic 
nuclear collisions at RHIC and LHC energies. These models are shown to predict qualitatively 
different associated hadron correlations with respect to tagged heavy quark jets. While both models 
feature similar far zone Mach and diffusion wakes, the far zone stress features are shown to be too 
weak to survive thermal broadening at hadron freeze-out. However, these models differ significantly 
in a near zone "Neck" region where strong chromo-fields sourced by the heavy quark jet couple to 
the polarizable plasma. Conical associated correlations, if any, are shown to be dominated by the 
jet induced transverse fiow in the Neck zone and unrelated to the weak far zone wakes. Unlike in 
AdS/CFT, we show that the induced transverse fiow in the Neck zone is too weak in pQCD to 
produce conical correlations after Cooper-Frye freeze-out. The observation of conical correlations 
violating Mach's law would favor the strongly-coupled AdS/CFT string drag dynamics, while their 
absence would favor weakly-coupled pQCD-based chromo-hydrodynamics. 

PACS numbers: 25.75.-q, 11.25.Tq, 13.87.-a 



I. INTRODUCTION 



Recent interest in Mach-like conical di-jet correlations 
[H is due to suggestions [1, Q that a measurement of the 
dependence of the cone angle on the supersonic jet ve- 
locity V could provide via Mach's law (cos 0m = Cs/v) a 
constraint on the average speed of sound in the strongly 
coupled Quark-Gluon Plasma (sQGP) [4] created at the 
Relativistic Heavy Ion Collider (RHIC). In Refs. |1| 
we explored the robustness of this interpretation using 
the string drag model of strongly coupled plasma-field 
interactions @, 0] in the context of the Anti-de Sit- 
ter/Conformal Field Theory (AdS/CFT) correspondence 
P . This AdS / CFT motivated model provides a detailed 
holographic description of the induced stress tensor Q 
in the wake of a heavy quark jet moving at a constant 
velocity through a static strongly-coupled A/" = 4 Su- 
persymmetric Yang-Mills (SYM) background plasma at 
finite temperature Tq. Direct tests of this AdS/CFT 
string drag model using the ratio of bottom to charm 
nuclear modification factors in high energy nuclear col- 
lisions at RHIC and LHC have been proposed in Ref. 
p^ . In this work, we concentrate on another observable: 
the hadron correlations associated with tagged identified 
heavy quark jets. 

The AdS /CFT stress solution for a supersonic heavy 
quark 0, [HI features the expected far zone Mach cone 
as well as a strong forward moving diffusion wake. The 
far zone response is well described in the strong coupling 
limit of an A/" = 4 SYM plasma by a "minimal" shear 
viscosity over entropy density ratio rj/s = l/47r (near 
the uncertainty principle limit p^). However, in [tI. 
it was noted that the strong forward diffusion wake in the 
far zone of the AdS/CFT solution could spoil the double- 



shoulder signature of the Mach wake in accord with the 
general discussion of far zone hydrodynamics in [3j. 

In it was shown that in the strict supergravity limit, 
Nc I, g'sYM ^ 1 but A = gsYM^c > 1, in fact the 
far zone wakes have such small amplitudes that they only 
lead to a single broad peak in the away-side hadronic 
correlation after Cooper-Frye (CF) freeze-out of the fluid 
p^ . However, the AdS/CFT string drag solution fea- 
tures a novel nonequilibrium "Neck" near zone, where 
especially strong transverse flow relative to the jet axis 
induces an apparent conical azimuthal correlation of as- 
sociated hadrons even after CF thermal broadening at 
freeze-out. 

In Ref. [H], the Neck zone was defined as the region 
near the heavy quark where the local Knudsen num- 
ber [3 is KNiX) = r, |V • M|/|M| > 1/3. where 
M'{X) = T°'{X) is the momentum flow field of mat- 
ter and Ts = 4?7/ (SsTq) > 1/ (SttTq) is the sound atten- 
uation length, which is bounded from below for ultra- 
relativistic systems [l^, [13 • In the Neck region the in- 
duced transverse flow is surpringly large and is unrelated 
to the far zone Mach's law. In AdS/CFT, this is the 
field-plasma coupling zone where the stress tensor has a 
characteristic interference form dependence on the coor- 
dinates, 0{^XTg/R'^) [13, m, with R denoting the dis- 
tance to the heavy quark in its rest frame. In contrast, 
the stress in the far zone has the characteristic 0(Tq ) 
form. In addition, very near the quark the self Coulomb 
field of the heavy quark contributes with a singular stress 
0(VA/i?4) 9]. 

The above strong coupling AdS/CFT results moti- 
vated us in the present work to study whether similar 
novel near zone field-plasma dynamical coupling effects 
arise in weakly coupled perturbative Quantum Chromo- 
dynamics (pQCD). In Refs. [l^, [13i[13| the heavy quark 
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jet induced stress in a weakly-coupled QGP (wQGP in 
contrast to sQGP) was computed analytically in the lin- 
ear response approximation based on the Asakawa-Bass- 
Miiller (ABM) |22|. [23t generalization of chromo-viscous 
hydrodynamics [24| . The ABM generalization concen- 
trates on the "anomalous diffusion" limit, where the con- 
ductivity is dominated by field rather than stochastic dis- 
sipative scattering dynamics. 

As in the AdS/CFT string drag model, the generic far 
zone Mach and diffusion wakes are also clearly predicted 
in the pQCD based ABM formulation [l^, [2l[ as we show 
for the case of a w = 0.9 heavy quark jet in Fig.l. How- 
ever, the question of whether the far zone Mach cone flow 
correlations survives CF freeze-out of the plasma was not 
addressed up to now. In this paper, we extend the work 
of [19,,^] by solving numerically the full nonlinear 3-l-lD 
relativistic hydro dynamic equations using the SHASTA 
hydro code [23|, supplemented with the chromo-viscous 
stress source derived in Refs. [19, 20,]. We specialize to 
the ideal fluid case of vanishing viscosity to minimize 
the dissipative broadening of any conical correlations and 
therefore maximizing the signal to noise ratio. 
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FIG. 1: (Color online) The fractional energy density pertur- 
bation Ae/eo = s{xi, Xp)/eo ~ 1 (in the lab frame) due to 
a heavy quark with v = 0.9 in a QCD plasma of tempera- 
ture To = 200 MeV. The induced fluid stress was calculated 
using 3-1- ID hydrodynamics [25j with the anomalous pQCD 
source of Neufeld (left: panel) and AdS/CFT ^ (right 
panel). A trigger jet (not shown) moves in the —x direc- 
tion. The away-side jet moves in x direction and contours 
of Ae/eo = —0.15,0.1,0.2,0.5,0.7 are labeled in a comoving 
coordinate system with xi — x — vt and the transverse radial 
coordinate Xp in units of l/nTo ~ 0.3 fm after a total transit 
time t — 5fm/c= 14.4/ (ttTq). The ideal Mach cone for a point 
source is indicated by the yellow dashed line in the xi — Xp 
plane. See Fig. [2] for a zoom of the Neck region inside of the 
black box. 



We emphasize that our aim here is not to address the 



current light quark/gluon jet RHIC correlation data. Our 
goal is to point out the significant differences between 
weakly coupled and strongly coupled models mechanisms 
of heavy quark energy loss that can be tested experi- 
mental when identified heavy quark (especially bottom 
quark) jet correlations will become feasible to measure. 
We limit this study to the most favorable idealized con- 
ditions (uniform static plasma coupled to the external 
Lorentz contracted color fields). Common distortion ef- 
fects due to evolution in finite expanding plasma geome- 
tries will be reported elsewhere. 

We use natural units below and Lorentz indices are 
denoted with Greek letters fj,,i> = 0, . . . , 3 while inter- 
nal SU{Nc) adjoint color indices are a — 1, . . . , ~ 1. 
Also, the Minkowski metric — diag(— , +) is 
employed. In our system of coordinates, the beam is in 
the z direction, the associated heavy quark jet moves 
along the x direction with velocity v = v x. We define 
xi — X — vt and Xp is the transverse cylindrical radial 
coordinate perpendicular to the jet axis. 
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FIG. 2; (Color online) A magnified view of the near "Neck" 
zone shows the relative local energy density perturbation 
Ae / eo and fluid flow directions induced by a heavy supersonic 
quark jet moving with v = 0.9. As in Fig.l, the pQCD con- 
tours were computed using 34- ID hydrodynamics 25] sourced 
by [131 (left panel). The AdS/CFT Neck zone 5] (right panel) 
uses numerical tables from 7] . The purple dashed line indi- 
cates the ideal far zone point source shock angle. The heavy 
quark is at the origin of these comoving coordinates. The 
arrows indicate both direction and relative magnitude of the 
fluid flow velocity. The numbers in the plot label the con- 
tours of constant Ae/eo. Note that Ae/eo is larger in pQCD 
but that the transverse flow generated near the quark is much 
stronger in the AdS/CFT model. 



II. STRESS ZONES 

Energetic back-to-back jets produced in the early 
stages of a heavy ion collision transverse to the beam 
axis traverse the sQGP and deposit energy and momen- 
tum along their path in a way that depends on the non- 
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equilibrium details of the physics of the field-plasma cou- 
pling. In the case when one of the jets is produced near 
the surface (the trigger jet), the other supersonic away- 
side jet moves through the plasma and generates in the 
far zone a Mach like conical perturbation in the plasma 
stress as seen in Fig.l. The resulting conical correlation 
(with respect to the away-side jet axis) is naturally ex- 
pected to lead to an enhancement of associated hadrons 
at the characteristic Mach angle jH, 01- Below some 
transverse momentum saturation scale [28|, jet physics 
depends on the properties of the medium and one can 
test different models of jet-medium coupling dynamics by 
studying the detailed angular and rapidity correlations. 
In this work we demonstrate a striking difference between 
strongly coupled AdS/CFT and moderate coupling, mul- 
tiple collision pQCD transport models for that coupling, 
and study their experimentally neasurable consequences. 

We ignore here the near-side associated correlations 
and focus on away-side jet-hadron azimuthal correlations. 
The energy-momentum stress induced by the away-side 
heavy quark jet in both pQCD and AdS/CFT can be con- 
veniently decomposed into four separate contributions as 
in 



(1) 

The static isotropic background stress tensor is assumed 
to be T^g = diag(e,p,p,p), where e = KTq is the corre- 
sponding background energy density of a gas of massless 
5t/(3) gluons Kqcd = 87rVl5 whereas for SU{Nc) SYM 
KsYM = 37r^(iV^ — l)/8. In both cases, e — 3p and the 
background temperature is Tq. 

The Coulomb contribution to the energy-momentum 
tensor ST^'^^^i{X) arises from the near zone Lorentz con- 
tracted Coulomb field that remains attached to the heavy 
quark since we consider only moderate but supersonic ve- 
locities Cs < V < 0.9. We can therefore neglect radiative 
energy loss that dominate in the ultrarelativistic case. 
The bare comoving Coulomb self-field stress has the sin- 
gular form ST^^'^^i cx in the quark rest frame. In 
both pQCD and AdS/CFT cases we subtract this vac- 
uum self-field stress as in Q. In other words, the zero 
temperature contribution to the in-medium stress ten- 
sor is always subtracted. While in AdS/CFT the form 
of the Coulomb tensor is known exactly @, in pQCD 
this contribution can only be calculated perturbatively. 
The leading order expression for the chromo-fields pro- 
duced by the source in pQCD, in the limit where the 
dielectric functions are set to unity, displays the same 
Lienard-Wiechert behavior as in AdS/CFT. 

The far zone "Mach" part of the stress can be expressed 
in terms of the local temperature T{X) and fluid flow ve- 
locity fields C/"(X) through the first-order Navier-Stokes 
stress form 



sT 



where s{X) is the local entropy density and d^'^U'^^ is 
the symmetrized traceless shear flow velocity gradient. 
In AdS/CFT, ri/s = l/47r [H. The theta function in Eq. 
defines the "far zone" that includes the Mach and Dif- 
fusion linearized hydrodynamic sound waves. In the far 
zone equilibration rate is over three times the local stress 
gradient scale, and first-order Navier-Stokes dissipative 
hydrodynamics provides an adequate description of the 
evolution in that zone. In AdS/CFT, the hydrodynamic 
description was shown to be valid down to distances of 
roughly 3/ (ttT) [16| (see also [sO] for an analysis of the 
far zone). The near Neck zone within a thermal Compton 
wavelength of the heavy quark is a generally nonequilib- 
rium region strongly influenced by the coupling of the 
heavy quark's bare classical Coulomb field to the polar- 
izable plasma. 



III. THE PQCD SOURCE TERM 

We use the nonlocal pQCD chromo-viscous stress 
source, S'^{x), derived in Ref. [l^l, to drive the per- 
fect fluid response of a pQCD fluid assuming the ideal 

= in order to maximize any pQCD transport induced 
azimuthal conical signature. Finite viscosity of course 
washes out some of the induced correlations as shown in 
PH . However, our main finding below is that even in 
this perfect 77 = hydrodynamic limit the induced corre- 
lations sourced in pQCD are too weak to generate conical 
correlations after freeze-out. 

Thus, we consider only the 77/s = limit of the 
full anomalous chromo-viscous equations derived in 
[23| . However, we retain the anomalous diffusion stress 
Neufeld source (Eqs. (23-24) of [20|). We can rewrite Eqs. 
(6.2 - 6.11) of [2^ in the more familiar covariant Joule 
heating form 



pua a ja 



{F^^'^cJcM^F^''") (3) 



where F^'^°(X) is the external Yang-Mills field tensor 
and r{X) = / d'^K/{2T:feyi]i{-iK ■ X)r{K) is the 
color current that is related via Ohm's law to F'^'^°'{K) 
through the (diagonal in color) conductivity rank three 
tensor a^ap{K). The * denotes a convolution over the 
nonlocal non-static conductive dynamical response of 
the polarizable plasma. The covariant generalization 
of Neufeld's source is most easily understood through 
its Fourier decomposition, J^iK) = a^fj,JK)F''°"'{K), 
with the color conductivity expressed as |3ll | 



Pf^Pa dp 



P-K + iP- U/t 



:MP) (4) 



rpfJ.V 



0(1 - ai^iv) 



where /o(P) = 2{N^-1)G{P) is the effec- 
tive plasma equilibrium distribution with G{P) — 
{2TT^)-^e{Po)6{P^)/{eP°^^ - 1). Here, t/^ is the 4- 
velocity of the plasma as in Eq. 12]). For an isotropic 
plasma U^af^apiK) = —<^^j.a{K)- In the long wave- 



(2) length hmit, U,}a'"^'^{K 0) 



^i)5'^"/3, where 
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mjj = g T is the Debye screening mass for a nonin- 
teracting plasma of massless SU{3) gluons in thermal 
equilibrium. 

The relaxation or decoherence time r* in (|4]) has the 
general form noted in [2^ r* = (1/rp + 1/tc + l/ran)^^ 
where Tp oc (a^T ln(l / as))~^ being the collisional mo- 
mentum relaxation time [HillSi ''c = {c( s N cT \n.{l / g))~^^ 
being the color diffusion time defined in (3l| . and Tan oc 
(to£i(77|V-U|/Ts)^/^)~^ being the anomalous strong elec- 
tric and magnetic field relaxation time derived in Eq. 
(6.42) of [2^. We note that one can express 
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1 



(5) 



in terms of the local Knudsen number = Tg/L used in 
Eq. (21) . Here L is the characteristic stress gradient scale. 
However, because ry cx t*sT, Eq. ([5]) is really an implicit 
equation for Tan- Combining these relations and taking 
into account the uncertainty principle constraint [l3| that 

bounds r* ^ 1/ (3T) for an ultrarelativistic (conformal) 
plasma, we have 



^ cx r(ai g'^lng ^ + a2g'^lng ^ + a^, g^jK^) ^ iT 

(6) 

where ai , a2 , 03 are numerical factors. In the near zone 
(see Fig. 2) close to the quark. Km gets large. Tan can 
become the dominant contribution to t* in the presence 
of strong classical field gradients. 

We would like to emphasize an important subtle point 
in the application of Eq. ^ to our heavy quark jet prob- 
lem. In order to neglect viscous dissipation in the pQCD 
response, the relaxation rate must be very large com- 
pared to the characteristic gradient scale. Hence, in the 
far zone at least the imaginary part of the conductiv- 
ity denominator in Eq.Q must be large and dominant. 
However, in the Neck region the field gradients become 
very large and the relevant wave numbers of the hydro 
response K 3> ST exceed the uncertainty limited equili- 
bration rate. Because we only need to consider the con- 
ductivity in the asymptotic large K limit in the near zone, 
it becomes possible to neglect the ~ i3T maximal relax- 
ation rate in the energy denominator and to formally set 
1 /r* — > 0^ - as if the coupling were parametrically small 
(as assumed in Eqs. (53-56) of 20]). Only in this high 
frequency, high wave number limit, relevant for the Neck 
zone physics, is the color conductivity computable as in 
[20| . The neglect of dissipation in the Neck zone maxi- 
mizes the acceleration of the plasma partons, which can 
subsequently generate transverse collective plasma flow 
relative to the jet axis. What we must check numeri- 
cally is whether this maximum transfer of field energy- 
momentum from the field to the plasma is sufficiently 
anisotropic to generate a conical correlation of the asso- 
ciated hadron fragments. 



IV. FREEZE-OUT PROCEDURES 

As noted previously, we consider here only the ideal- 
ized static medium here to maximize plasma response 
signals. Distortion effects due transverse expansion and 
non-conformal equation of states with a phase transi- 
tion will be presented elsewhere. These effects, while 
important for phenomenological comparisons to heavy 
ion data, however obscure the fundamental differences 
between weak and strong coupled dynamics that is our 
focus here. 

Given the large theoretical systematic uncertainty in- 
herent in any phenomenological model of nonperturba- 
tive hadronization, we consider here two simple limits 
for modeling the fluid decoupling and freeze-out. In one 
often used limit, we freeze-out computational fluid cell 
via the CF prescription on an isochronous hypersurface. 
This scheme takes into account maximal thermal broad- 
ening effects. In the opposite limit, we freeze-out as- 
suming isochronous sudden breakup or shattering of fluid 
cells conserving only energy and momentum and avoid- 
ing hadronization altogether as described in more detail 
below. The difference between the two schemes provides 
a measure of the systematic theoretical uncertainty asso- 
ciated with the unsolved problem of hadronization. 

In the CF method, the conversion of the fluid into free 
particles is achieved instantaneously at a critical sur- 
face dYjg lldl . If we assume such a freeze-out scheme 
[li S Ha l32| |. the particle distributions and correlations 
can be obtained from the flow velocity field U^{X) and 
temperature profile T{X). For associated (massless) par- 
ticles with — {pt,pt cos{tt — (/)) , pr sin^n — (/)),0) the 
momentum distribution at mid rapidity y — is 



dN 



PTdpTdydcf) 



= [ dI]^pq/o(^7^p^T)-/e,] (7) 

V=0 Js 



where px is the transverse momentum, S(X) is the 
freeze-out hypersurface, and /o — exp{—U^Pi^/T{X)) is 
a local Boltzmann equilibrium distribution. No viscous 
corrections to Eq. ^ are included since we are working 
here in the perfect fiuid limit with i] = 0. We subtract the 
isotropic backgroundyield via feq = f\ui^=o,T=To- More- 
over, we follow 0, [1[ [2^ and perform an isochronous 
freeze-out where dS^ = x (1, 0, 0, 0). 

We remark that the absence of well-defined quasi- 
particle states in AdS/CFT plasmas at large t'Hooft cou- 
pling indicates that CF can only, at best, give a quali- 
tative idea of the observable hadron level angular cor- 
relations 0- Moreover, even in the pQCD quasiparti- 
cle limit, CF freeze-out remains a strong model assump- 
tion. In the pQCD case, in the associated pr range of 
interest a coalescence/recombination hadronization sce- 
nario [ssl . [34] may be more appropriate. However, we 
expect similar CF thermal broadening effects if coales- 
cence hadronization is assumed and full three momentum 
conservation is taken into account. 

As an alternate freeze-out scheme we consider a 
calorimetric-like observable given by the momentum den- 
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sity weighted polar angle distribution relative to the jet 
axis: 
dS 



dcost 



|^c|'5(cos6l- cos 



cells 



/.»x|M(X)|.(».«-^) 



(8) 



This quantity differs from CF mainly by the neglect of the 
thermal smearing at the freezeout time, and thus it max- 
imally amplifies the angular anisotropy of the associated 
hadrons. The very strong assumption in this decoupling 
scheme is that hadrons from each frozen-out cell emerge 
parallel to the cell total momentum = d'^x T''*(x, tf). 
Here 9 = tt — Otrigger is the polar angle with respect 
to the away-side heavy quark jet. Many other similar 
purely hydrodynamic measures of bulk flow are possible 
[3], e.g. entropy instead of momentum density weighted. 
However, we found no qualitative differences when the 
weight function is changed. We used a narrow Gaus- 
sian approximation to the Dirac delta in Eq. ([8]) with a 
A cos 6 — 0.05 width and checked that the results did not 
change significantly if the width was varied by 50%. 

Our goal in this letter is not to decide which hadroniza- 
tion scheme is preferred but to apply these two com- 
monly used measures of hydrodynamic to help quan- 
tify the observable differences between two strongly dif- 
ferent approaches to jet-plasma interactions (pQCD x 
AdS/CFT). The CF freeze-out employed here and in 
Refs. d, [1^, [131 is especially questionable in the non- 
equilibrium Neck region but provides a rough estimate of 
intrinsic thermal smearing about the local hydrodynamic 
flow. In Fig. 2 we show the relative local energy distur- 
bance and flow proflle in the Neck region created by a 
v=0.9 jet in both pQCD and AdS/CFT. Note that the 
relative transverse flow in the Neck zone in AdS/CFT is 
significantly larger than in pQCD and as we show below 
this is reflected in the final angular correlations from that 
region in both hadronization schemes. 



V. FREEZE-OUT RESULTS IN PQCD 

The initial away-side heavy quark jet is assumed to 
start t — aX xi = —4.5 fm and the freeze-out is done 
when the heavy quark reaches the origin of the coordi- 
nates at time tf — 4.5/w fm. This provides a rough de- 
scription of the case in which a uniformly moving heavy 
quark punches through the medium after passing through 
4.5 fm of plasma. 

The numerical output of the SHASTA code is the tem- 
perature and fluid flow velocity flelds T{X) and \J{X). 
The hydrodynamic equations were solved in the pres- 
ence of the source term iS''(X) computed analytically 
by Neufeld in [2^ in the limit where the dielectric func- 
tions that describe the medium's response to the color 
fields created by the heavy quark were set to unity. The 
effects from medium screening on S'^ were studied in de- 
tail in Ref. i20l. In our numerical calculations we used 



Xpmax = I/tid as an infrared cutoff while the minimum 
lattice spacing naturally provided an ultraviolet cutoff. 
The background temperature was set to Tq = 0.2 GcV. 
We assumed Us — f/^/ (47r) — I/tt in our calculations 
involving the pQCD source. 

The results for the bulk flow according to Eq. ^ in 
pQCD are shown in the upper panel in Fig. [3) The curves 
are normalized in such a way that the largest contribu- 
tions are set to unity. Note that the pQCD bulk energy 
flow distribution has a large forward moving component 
in the direction of the jet for all the velocities studied 
here. In the far zone, this forward moving energy flow 
corresponds to the diffusion wake studied in . The red 
curve with triangles in the upper panel in Fig. [3] cor- 
responds to the yield solely from the Neck region for 
V = 0.9. The very small dip at small 6* = is mostly 
due to the weak Neck zone pQCD but the most of the 
momentum flow from both Neck and Diffusion wakes is 
directed along the jet axis. The relatively small trans- 
verse energy flow in the Neck region is evident on the left 
panel of see Fig. 2 in contrast to the much larger trans- 
verse flow predicted via AdS in that near zone. The Mach 
cone emphasized in Ref. [l^ [2l| is also clearly seen but 
its amplitude relative to the mostly forward diffusion plus 
Neck contribution is much smaller than in the AdS / CFT 
case. The weak Mach peak roughly follows Mach's law 
as V approaches c^. 

In Fig. 4, our CF freeze-out results for the associated 
away-side aziniuthal distribution for v — 0.58, 0.75, 0.9 
for mid-rapidity and pr = 57rTo ~ 3.14 GeV light 
hadrons are shown. The pQCD case , computed using 
the output of the SHASTA code into Eq. ([7|) , are shown 
in the upper panel. We show the angular function 
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dNic/,) 



Nmax \PTdpTdyd(f> 



y=o 



(9) 



where N^ax is a constant used to normalize the plots 
(note that this function is not positivc-deflnite). The 
pQCD angular distribution shows only a sharp peak at 
<f> = TT for all velocities. The red curve with triangles 
denotes the contribution from the pQCD Neck region for 
V = 0.9. Note that the different peaks found in the bulk 
flow analysis of the pQCD data shown in the upper panel 
in Fig. [3] do not survive CF freeze-out. We checked that 
no other structures appear if we either double pT to 5 
GeV or increase as to 0.5. We conclude that the strong 
forward moving diffusion wake as well as the mostly for- 
ward bow shock Neck zone dominate the away-side peak 
and that the thermal broadened Mach correlations are 
too weak in pQCD to contribute to the flnal angular cor- 
relations. 



VI. FREEZE-OUT RESULTS IN ADS/CFT 

We used the same setup employed in Q to perform 
the CF freeze-out of the A/" = 4 SYM AdS/CFT data 
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FIG. 3: (Color online) The (normalized) momentum weighted 
bulk flow angular distribution as a function of polar angle with 
respect to the away-side jet is shown for v = 0.58 (black), 
V = 0.75 (magenta), and v = 0.90 (blue) comparing pQCD 
anomalous chromo-hydrodynamics to the AdS/CFT string 
drag @, 01 model analyzed in Ref. [^. The red line with tri- 
angles represents the Neck contribution for a jet with v = 0.9 
and the arrows indicate the location of the ideal Mach-cone 
angle given by cos 6m = Cs/v, where Cs = 1/VS. 



computed by Gubser, Pufu, and Yarom in Ref. [3]. They 
calculated the energy-momentum disturbances caused by 
the heavy-quark, which in this steady-state solution was 
created at i — s- —oo and has been moving through the in- 
finitely extended A/" = 4 SYM static background plasma 
since then. The freeze-out is computed when the heavy 
quark reached the origin of the coordinates. The mass 
of the heavy quark M in the AdS/CFT calculations is 
such that M/Tq ^ -s/A, which allows us to neglect the 
fluctuations of the string [s^ . At Nc = 3 the simpli- 
fications due to the supergravity approximation are not 
strictly valid but it is of interest to extrapolate the nu- 
merical solutions to study its phenomenological applica- 
tions. We chose a plasma volume to be the forward light- 
cone that begins at xi = —4.5 fm and has a transverse 
size of Xp < 4.5 fm at Tq = 0.2 GeV (our background 
subtracted results do not change when larger volumes 
were used). Note that we assumed the same background 



temperature for both pQCD and AdS/CFT. The map- 
ping between the physical quantities in A/" = 4 SYM and 
QCD is a highly non trivial open problem (see, for in- 
stance, the discussion in Ref. 36]). We therefore again 
use CF and bulk momentum flow as two extreme limits 
to gauge possible systematic uncertainties. 

The (normalized) bulk momentum flow associated with 
the AdS/CFT data, computed using Eq. ([5]), is shown 
in the lower panel of Fig. [3l The bottom panel in Fig. 
[3] simply shows that in AdS/CFT there are more cells 
pointing in a direction near the Mach cone angle than 
in the forward direction (diffusion wake) when v = 0.9 
and V = 0.75 unlike in the pQCD case shown in the 
upper panel. However, when v = 0.58 the finite angle 
from the Mach cone is overwhelmed by the strong bow 
shock formed in front of the quark, which itself leads to 
small conical dip not at the ideal Mach angle (black ar- 
row). The red line with triangles in the bottom panel of 
Fig. [3] shows that the relative magnitude of the contribu- 
tion from the Neck region to the final bulk flow result in 
AdS/CFT is much smaller than in pQCD. However, note 
that small amplitude peak in the AdS/CFT Neck curve 
is located at a much larger angle than the correspond- 
ing peak in the pQCD Neck, as one would expect from 
the transverse flow shown Fig. 2. Moreover, one can see 
that a peak in the direction of the trigger particle can be 
found for all the velocities studied here. This peak rep- 
resents the backward flow that is always present vortex- 
like structures created by the jet as discussed in detail in 
Ref. [131 • Possible phenomenological consequences due to 
these vortex structures on the polarization of hyperons 
in heavy ion collisions were discussed in Ref. [13] ■ 

Our results for the CF freeze-out of the AdS/CFT solu- 
tion for V = 0.58, 0.75, 0.9 at mid-rapidity andpT — SttTq 
GeV in the lower panel in Fig. 21 A double peak structure 
can be seen for v = 0.9 and v — 0.75. Note, however, that 
the peaks in the AdS/CFT correlation functions do not 
obey Mach's law. This is because these correlations come 
from the Neck region where there is a strong transversal 
non-Mach flow [5|. This is explicitly shown by the red 
curve with triangles that represents the Neck contribu- 
tion for a jet with v = 0.9 as in Fig. [31 For v = 0.58, 
the resulting flow is not strong enough to lead to non- 
trivial angular correlations. A detailed study of the pr 
dependence of the away-side correlations associated with 
AdS/CFT heavy quark jets will be presented in a future 
work. In addition, the negative yield present in the CF 
curves for v — 0.58 and v — 0.75 is due to the presence 
of the vortices discussed above. 

In general [5|], the weak sound waves produced by a 
jet do not lead to a cone-like signal independently of the 
detailed flow and interference patterns because thermal 
smearing washes out the signal. Formally, if linearized 
hydrodynamics applies and in the low momentum limit 
(U • p << T), the associated hadron away-side distri- 
bution is only a very broad peak about <f> — n regard- 
less of the detailed combination of Mach wakes, diffusion 
wakes, or vortex circulation [1, [^, [23|. This result in- 
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FIG. 4: (Color online) Normalized (and background sub- 
tracted) azimuthal away-side jet associated correlation after 
Cooper-Frye freeze-out CF{(f>) (see Eq. ^ for pQCD (top) 
and AdS/CFT from (bottom). Here CF((^) is evaluated 
at pt ~ 5ttTo ~ 3.14 GeV and y — 0. The black line is for 
V = 0.58, the magenta line for v = 0.75, and for the blue 
line V = 0.9. The red line with triangles represents the Neck 
contribution for a jet with v = 0.9. 



volving CF freeze-out can only be circumvented either in 
regions with high flow velocities and large gradients as in 
the Neck zone [5[ , or by increasing px to unrealistic high 
values d,!!!]. 

One of the main differences between the two freeze-out 
procedures we employed (in both AdS/CFT and pQCD) 
concerns the relative magnitude of the contribution from 
the Neck region to the final angular correlations: the 
Neck region is much more important in CF than in the 
bulk flow measure computed via Eq. ([8]). This is due to 
the exponential factor in CF, which largely amplifies the 
contribution from the small region close to the jet where 
the disturbances caused by the heavy quark become rel- 
evant. 

VII. CONCLUSIONS 

In this paper, we showed that the angular correlations 
obtained after an isochronous Cooper-Frye freeze-out of 



the wake induced by punch-through heavy quark jets (in 
a static medium) in the Neufeld et al. pQCD model 
of anomalous chromo- viscous hydrodynamics do not dis- 
play a conical structure. This should be compared to 
the conical-like structures seen after CF freeze-out of the 
strongly-coupled AdS/CFT string drag model. We ex- 
pect similar thermal broadening effects to occur if alter- 
native coalescence/recombination hadronization models 
[S^ . [S^ l are used. The isochronous hypersurface we used 
is needed in order to compare AdS/CFT to pQCD since 
AdS/CFT heavy quark solutions have only been com- 
puted so far in a static medium. For realistic simulations 
that can be compared to data, effects from the medium's 
longitudinal, transverse, and elliptic flow must be taken 
into account as discussed in detail in [s^ . 

Unlike AdS/CFT, the conical flow from the associated 
nonequilibrium Neck zone in pQCD (see the red region in 
the left panel of Fig. 2 and the red curve in Fig. 3) is too 
weak to survive CF freeze-out. In both cases, the actual 
Mach wakes do not appear after standard CF freeze-out. 
Mach-like peaks are only observable in the sudden shat- 
tering freeze-out scenario described in Eq. ([5]) in both 
pQCD and AdS/CFT, in which thermal broadening is 
entirely neglected. 

The Neck region (in both pQCD and AdS/CFT) gives 
the largest contribution to the total yield in CF freeze-out 
while its contribution in the other extreme case involving 
the bulk flow hadronization is not as relevant. This in- 
dicates that the magnitude of the Neck's contribution to 
the flnal angular correlations is still strongly model de- 
pendent. Nevertheless, our results suggest that conical 
but non-Mach law correlations arc much more likely to 
appear in AdS/CFT than in pQCD. 

We propose that the measurement of the jet veloc- 
ity dependence of the associated away-side correlations 
with identified heavy quark triggers at RHIC and LHC 
will provide important constraints on possible pQCD 
versus AdS/CFT dynamical non-Abelian field - plasma 
(chromo-viscous) coupling models. 
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